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Abstract:  The  cohesion  of  a  grain  boundary  (GB)  is  believed  to  be  the  control¬ 
ling  factor  limiting  the  ductility  of  high-strength  metallic  alloys,  and  particularly 
those  containing  W.  Intergranular  embrittlement  is  usually  associated  with 
segregation  of  impurities  at  the  GBs.  Impurities  present  in  ppm  concentrations 
can  result  in  a  dramatic  decrease  in  plasticity.  This  paper  reviews  recent  results 
on  both  semi-empirical  and  first-principles  modelling  of  the  energetics  and  the 
electronic  structures  of  impurities  on  a  23  (111)  GB  in  W.  Our  calculations 
have  shown  that  impurities,  such  as  N,  O,  P,  S,  and  Si,  weaken  the  intergranular 
cohesion  resulting  in  ‘loosening’  of  the  GB.  The  presence  of  B  and  C  on  the 
contrary,  enhances  the  interatomic  interaction  across  the  GB.  The  so-called 
‘site-competition  effect'  should  play  an  important  role  affecting  impurity 
distribution  in  W  GBs.  Among  the  impurities  analyzed,  B  in  the  GB  has  the 
lowest  energy  and  thus  would  tend  to  displace  other  impurity  atoms  from  the 
GB.  Microalloying  with  10-50  ppm  B  may  be  an  effective  way  of  improving 
tungsten’s  ductility.  These  results  are  important  for  understanding  the 
fundamental  physics  of  intergranular  embrittlement. 


INTRODUCTION 

The  reduced  cohesion  of  grain  boundaries  (GBs) 
is  often  the  controlling  factor  limiting  ductility, 
and  hence  performance  and  reliability  of  high- 
strength  metallic  alloys.'  Intergranular  embrittle¬ 
ment  in  metals  is  usually  caused  by  impurities 
segregating  towards  the  GBs.^''’"  A  ductile-brittle 
transition  temperature  (DBTT)  as  low  as 
- 1 96°C’  was  observed  in  high-purity  W  single 
crystals  obtained  by  electron-beam  zone-melting 
with  impurity  gettering.  Impurities  present  in  bulk 
concentrations  of  only  10'^- 10"“*  atomic  percent 
can  result  in  a  dramatic  decrease  in  plasticity.  This 
drastic  degradation  of  mechanical  properties  of 
metallic  alloys,  in  particular  W,  poses  significant 
technological  and  application  problems.  This 
detrimental  effect  of  minute  impurity  concentra¬ 
tions  can  be  readily  understood;  a  ppm  amount  of 
impurity  is  sufficient  for  saturating  all  the  GBs  in  a 
polycrystal  of  a  typical  grain  size.  Sensitivity  of  the 
DBTT  to  the  grain  size  confirms  the  above  physi¬ 
cal  concept:  the  larger  the  grain  size,  the  smaller  is 
the  amount  of  impurity  needed  to  saturate  the 


GB.^  Finer-grain  polycrystals  are  known  to  be  less 
brittle.  For  a  23  (111)  GB  in  a  W  BCC  crystal. 
Table  1  shows  the  estimated  amounts  of  im¬ 
purities  sufficient  to  create  a  monolayer  coating. 

Brittle  fracture  is  inherently  related  to  the 
nucleation  and  propagation  of  cracks.  In  this 
paper,  however,  we  shall  focus  exclusively  on  the 
microscopic  aspects  of  embrittlement;  the  role  of 
impurities  will  therefore  be  considered  only  from 
the  point  of  view  of  their  effect  on  intergranular 
cohesion  as  the  main  factor  controlling  ductile/ 
brittle  behavior. 

If  the  presence  of  impurities  is  so  detrimental, 
gettering  these  impurities  is  an  obvious  way  of 
ductilizing  W.  A  well-known,  though  extremely 
costly,  option  is  to  use  the  so-called  ‘Rhenium 
Effect’  (see,  for  example.  Ref.  8).  A  more  promis¬ 
ing  way  of  removing  the  ‘harmful’  impurities,  such 

Table  1.  Amounts  of  impurity  atoms  sufficient  for  forming  a 
monolayer  in  the  GB 

Grain  size  (^m)  5  10  20  50  100 

Amount  of  impurity  (atomic  ppm)  60  30  15  6  3 
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as  O,  N,  P,  from  the  GBs  is  by  gettering  where 
these  form  thermodynamically  stable  phases  with 
other  elements,  e.g.  Ti,  Y,  Mo,  Zr,  Hf,  This 
gettering  approach,  however,  requires  precise 
control,  since  any  improved  ductility  will  occur 
only  as  long  as  these  second-phase  precipitates 
remain  fine;  any  excessive  growth  of  precipitates, 
such  as  the  so-called  ‘Ostwald  ripening’,  would 
result  in  an  adverse  embrittling  effect  due  to  the 
formation  of  these  latent  crack  starters. 

During  the  recent  decade,  extensive  experi¬ 
mental  work  was  directed  at  a  better  understand¬ 
ing  of  the  effect  of  impurities  on  brittleness  of  W. 
In  this  respect,  a  significant  contribution  of  Rus¬ 
sian  metallurgists  should  be  acknowledged  (see, 
for  example,  refs  7,  9-14  and  references  therein). 
Unfortunately,  most  of  the  related  papers  have 
been  published  in  Russian  and  are  therefore  virtu¬ 
ally  unknown  to  metallurgists  in  the  West,  though 
some  of  the  papers  have  been  translated. 

Recent  progress  in  developing  efficient 
methods  of  first-principles  calculations  and  com¬ 
putational  algorithms  made  possible  systematic 
studies  of  the  role  of  impurities  in  intergranular 
cohesion  of  transition  metals  at  the  atomic  and  the 
electron-ion  level.  Calculations  using  cluster,  two- 
dimensional,  and  supercell  models  of  GBs  with 
impurities  have  provided  an  in-depth  insight  into 
mechanisms  of  GB  decohesion  (for  references, 
see  ref.  15). 

Since  the  first-principles  electronic  calculations 
on  low-symmetry  systems  (such  as  lattice  defects 
or  GBs)  are  still  extremely  complicated  and  costly, 
semi-empirical  methods  based  on  strong  first- 
principles  foundations  have  also  been  developed. 
Among  them,  the  most  popular  is  the  Embedded 
Atom  Method  (EAM).*'’  This  method  has  been 
successfully  used  in  a  wide  variety  of  calculations. 

The  purpose  of  this  paper  is  to  elucidate,  at  the 
electron-atom  level,  the  energetics  of  impurities 
in  a  tungsten  GB  and  analyze  the  effect  of  im¬ 
purities  on  W  intergranular  cohesion.  A  deeper 
understanding  of  the  cohesion-decohesion  pro¬ 
cesses  at  the  microscopic  level  will  lay  a  founda¬ 
tion  for  ’smart  design’  of  ductile  W  alloys.  In 
particular,  the  theoretical  analysis  of  the  electron 
structure  and  the  energetics  of  W  GBs,  both  clean 
(CL)  and  with  impurities,  enables  one  to  make 
important  predictions.  As  a  result  of  this  theore¬ 
tical  analysis,  we  suggest  a  way  of  improving  the  W 
ductility  by  using  the  so-called  ‘site-competition’ 
effect.  Boron  introduced  in  minute  quantities  of 
10-50  atomic  ppm  could  cleanse  the  W  GBs  of 
other  harmful  impurities,  enhancing  the  inter¬ 


granular  cohesion  and  thus  improving  the  ductil¬ 
ity. 


ENERGETICS  OF  IMPURITIES  IN  W  GBs 

AND  THE  SITE-COMPETITION  EFFECT 

In  order  to  study  the  energetics  of  impurity  atoms 
in  a  W  GB,  we  have  chosen  first  to  calculate  the 
quantity  that  may  be  called  the  ‘environment- 
sensitive  embedding  energy’  (ESE):  the  energy  of 
an  impurity  atom  in  an  atomic  environment 
typical  for  a  GB.  Knowledge  of  these  energies  for 
various  impurities  enables  one  to  compare  the  .i 

relative  stability  of  a  particular  impurity  in  the 
W-GB  environment. 

Having  calculated  the  ESEs  for  a  number  of  , 

impurity  atoms,  one  can  use  this  information  in  a 
modified  EAM  approach  for  calculating  the  GB 
relaxation.  The  latter  calculation  enables  one  to 
draw  important  conclusions  regarding  the  inter¬ 
granular  cohesion  in  W  in  the  presence  of  a 
specific  impurity  in  the  GB. 

The  model  chosen  for  the  GB  environment  is  an 
eight-atom  hexagonal  supercell  (W^,#,  where  •  is 
an  impurity  atom).  The  supercell  is  shown  in  Fig. 

1,  together  with  the  capped  trigonal-prism  co¬ 
ordination  of  the  surrounding  W  atoms.  A  tri¬ 
gonal-prism  GB  configuration  is  believed  to  be  a 
typical  GB  environment  in  BCC  metals  and  is 
predicted  by  the  theory  of  hard-sphere  packing. 
Atomistic-relaxation  studies  have  shown  that,  in 
BCC  Fe,  an  impurity  atom,  such  as  P  or  B,  is  likely 
to  occupy  an  interstitial  position  in  the  center  of 
the  trigonal  prism  formed  by  Fe  atoms  in  the  GB 
core  (even  if,  as  in  the  case  of  P  and  B,  the  im¬ 
purity  forms  a  substitutional  solid  solution  with 
the  host).  The  hexagonal  superceU  of  Fig.  1  has  a 
relatively  high  symmetry;  it  also  emulates  a  ( 1 1 1 ) 

23  GB  environment.’^ 

We  performed  the  spin-polarized  scalar-relati-  • 

vistic  Linear  Muffin  Tin  Orbitals  (LMTO)  cal¬ 
culations  (our  method  and  approximations  were 
the  same  as  in  ref.  15).  First,  a  series  of  calcula-  * 

tions  (for  six  different  volumes)  was  performed 
with  an  impurity  absent  from  the  supercell,  i.e.  an 
empty  sphere  of  the  same  radius  as  that  of  the 
impurity’s  atomic  sphere  was  substituted  for  the 
latter.  Similar  calculations  were  then  performed 
for  each  of  the  impurities:  B,  C,  N,  O,  Al,  Si,  P, 
and  S.  The  ESEs  were  defined  as  follows: 

ESE  =  E( W,#)  -  £( W,0)  -  £(•)  ( 1 ) 

where  E(W^,#)  and  E(Wf,0)  are,  respectively,  the 
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Fig.  1.  The  W^#  hexagonal  supercell  emulating  a  typical  trigonal- prism  environment  of  W  atoms  in  the  ( 1 1 1 )  23  FB;  (a)  the 
supercell;  (b)  the  trigonal-prism  co-ordination;  ®  W,  •  Impurity.  Numbers  designate  atoms  W1  to  W4. 


energies  of  the  supercell  with  and  without  the 
impurity  (O  stands  for  an  empty  sphere  sub¬ 
stituted  for  the  impurity  atom),  and  E(#)  is  the 
energy  of  the  free  impurity  atom.  In  order  to  make 
the  calculations  more  consistent,  we  have  chosen 
to  use,  as  E(#)s,  the  values  of  E(Wft#)-E(W^O) 
extrapolated  to  the  zero-charge  density  (n  =  0), 
which  would  correspond  to  the  energy  of  an 
impurity  in  the  GB  environment  with  the  host 
crystal  lattice  infinitely  expanded.  The  ESE 
energies,  eqn  ( 1 ),  as  a  function  of  n,  the  electron- 
charge  density  due  to  W  atoms  at  the  impurity 
site,  are  presented  in  Fig.  2. 

Plots  in  Fig.  2  explain  an  experimentally 
observed  phenomenon  known  as  the  ‘site-com¬ 
petition'  effect.  As  one  can  see,  in  the  range  of 
electron-charge  density  typical  of  a  GB 
(0-01 5-0-025  a.u.),  B  has  the  lowest  energy  and 
thus  would  tend  to  displace  the  other  impurities 
off  the  GB.  Thus  there  exists  a  ‘site-competition 
hierarchy’.  In  fact,  in  W,  N  was  found  to  compete 
successfully  over  C,' '  while  C  ‘beats’  P.’’ 

The  plots  in  Fig.  2  also  reveal  an  important 
aspect  of  GB-impurity  behavior.  All  the  plots 
have  well-pronounced  minima.  The  positions  of 
the  minima  corresond  to  the  electron  density  at 
the  impurity  site  due  to  the  surrounding  W  atoms 
that  would  occur  if  the  GB  were  allowed  to  relax 
in  such  a  way  as  to  minimize  the  impurity’s 


0.00  0.01  0.02  0.03 

n(a.u.) 


Fig.  2.  The  ‘environment-sensitive  embedding  energies' 
plotted  against  n,  the  electron-charge  density  {in  atomic 
units,  a.u.). 

energy.  The  minima  positions  systematically 
(except  for  N)  shift  towards  lower  densities  as  the 
impurity  loses  its  competitive  power.  A  smaller 
charge  density  means  a  ‘looser’  GB,  less  strong 
and  more  prone  to  decohesion.  The  minimization 
of  the  total  GB  energy  (rather  than  only  the 
energy  of  the  impurity  atom)  gives  the  charac¬ 
teristic  charge  densities,  which  are  somewhat 
higher  than  those  at  the  minima.  However,  from 
this  point  of  view,  O,  S,  P,  Si,  and  AI  are  the 
obvious  ‘decohesion’  candidates,  whereas  B,  C, 
and  N  may  be  labeled  ‘cohesion  enhancers’.  In 
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fact,  B  and  C  were  experimentally  found  to 
improve  the  GB  cohesion  in  whereas 

O,  Si,  P,  and  S,  known  as  strong  embrit- 
tlers,’-'’  '^'-'  are  believed  to  weaken  the  GB  cohe¬ 
sion. 

As  mentioned  above,  the  GB  environment  we 
were  dealing  with  was  that  of  the  ( 1 1 1 )  23  tilt  GB. 
The  GB  structure  can  be  represented  as  a  succes¬ 
sion  of  ( 1 1 1 )  hexagonal  planes: 

. . .  CBACBACB  ACB  ABC  ABC  ABC  ABC . . . 

(the  GB  plane  is  marked  by  A).  The  CBABC 
atomic  stacking  of  the  core  of  the  GB  (CL  or  with 
an  impurity)  is  just  the  one  emulated  by  the  eight- 
atom  supercell  shown  in  Fig.  1.  In  order  to  find 
the  GB  structure  corresponding  to  a  minimum  of 
energy,  eqn  (2),  the  interplanar  distances  were 
varied,  whereas  the  interatomic  spacings  and  the 
structure  within  the  (111)  planes  were  left 
unchanged. 

The  total  energy,  E,  was  calculated  by  using  a 
modified  EAM  approach: 

£  =  2«£emb(n(i?))+(l/2)2«,«-F( /?,/?') 

+  ESE(n(/?;„,p)) 

where  and  V{R,R')  -axe  the  EAM 

embedding  energy  and  the  pair  potential  as  found 
for  the  bulk  BCC  W  (we  used  the  Finnis-Sinclair 
functions  and  parameters).--  The  third  term  is  the 
energy  of  the  impurity  atom.  The  parameters  R 
and  R'  are  the  positions  of  the  host  atoms,  R^^^  is 
that  of  the  impurity,  and  n{R)  and  n{R-^^p)  are  the 
electron-charge  densities  at  the  site  of  a  host  atom 
and  the  impurity,  respectively.  The  electron- 
charge  density  at  a  given  site  was  taken  to  be  a 
superposition  of  the  free-atom-charge  densities. 


As  follows  from  our  calculations  (Table  2),  the 
volume  difference  between  the  CL  GB  and  that 
with  an  impurity  increases  monotonically  from  B 
through  S.  From  the  intuitive  point  of  view,  the 
increase  in  volume  is  expected  to  result  in  a 
weakening  of  interatomic  bonds,  though,  in  prin¬ 
ciple,  an  impurity  may  exert  a  stronger  interaction 
in  spite  of  the  lattice  expansion. 

Table  2  also  shows  the  interatomic  distances  in 
relaxed  GBs.  The  atom  designation  is  shown  in 
Figs  1,  3,  and  4.  We  also  number  the  atomic 
planes  (see  below)  accordingly:  plane  1  (the  GB 
plane)  contains  atoms  Wl,  plane  2  atoms  W2,  etc. 
The  distance  between  an  impurity  atom  (or  an 
empty  sphere,  O)  and  Wl  atom  (2-584  A),  is  the 
same  in  all  cases,  since  atoms  within  the  (111) 
planes  were  not  allowed  to  relax. 

The  counterpart  of  the  nearest-neighbor  dis¬ 
tance  in  the  bulk  BCC  lattice  (2-741  A)  is  the  dis¬ 
tance  between  two  W  atoms  in  the  ( 1 1 1 )  direction 
(W1-W4).  In  the  CL  GB,  this  distance  is  larger 
than  in  the  bulk,  while  the  shortest  distance  is  the 
one  between  two  W2  atoms  (W2-W2)  across  the 
GB.  In  the  CL  GB,  there  is  a  significant  void 
(occupied  by  an  O);  the  distance  between  W3 


Fig.  3.  Arrangement  of  schematic  atomic  planes  in  GB 
cores:  (a)  CL  GB;  (b)  GB  with  an  impurity  atom. 


Table  2.  Interatomic  distances  (in  A)  in  relaxed  GBs,  and  GB  volume-relaxation  effect  (A')* 


# 

Distances 

Volumes 

W1-W4 

WI-W2 

W2-W2 

Q 

Q-Q„ 

0 

1-807 

2-841 

2-873 

2-512 

148-63 

5-95 

0 

B 

2-170 

3-056 

2-914 

2-690 

155-28 

12-60 

6-65 

C 

2-223 

3-090 

2-924 

2-736 

156-48 

13-80 

7-85 

N 

2-290 

3-133 

2-940 

2-805 

158-12 

15-44 

9-49 

O 

2-317 

3-150 

2-948 

2-839 

158-80 

16-12 

10-17 

A1 

2-539 

3-303 

3-035 

3-182 

165-10 

22-34 

16-47 

Si 

2-545 

3-308 

3-037 

3-190 

165-27 

22-41 

16-64 

P 

2-588 

3-345 

3-048 

3-234 

166-37 

23-69 

20-44 

S 

2-666 

3-426 

3-090 

3-389 

169-07 

26-39 

20-44 

is  the  volume  of  nine-layer  GB,  is  that  of  the  CL  GB,  and  Q^cc  is  the  bulk  volume  per  nine  atoms,  142-86  A^  •  is  an 
impurity  atom,  O  is  an  empty  sphere  (vacancy). 
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group  # 

Fig.  4.  A£  =  £gb(*)“£gb(CL),  the  energy  difference 
between  the  GB  with  impurity  •  and  CL  GB,  plotted  against 
the  periodic-chart  group  number. 

atoms  (W3-W3)  across  the  void  (and  the  GB 
plane)  is  quite  large;  3-614  A.  Thus,  in  the  CL 
GB,  the  strongest  interaction  is  W2-W2,  followed 
by  Wl-Wl  and  W1-W2.  With  an  impurity  atom 
in  place  of  O,  the  interatomic  interaction  changes 
significantly.  Now,  the  shortest  distance  is  that 
between  an  impurity  atom  and  W3,  •-W3,  the 
distance  between  W2  atoms  across  the  GB  being 
larger.  The  interaction  between  the  impurity  atom 
and  atom  W3  becomes  of  utmost  importance.  It  is 
actually  responsible  for  the  intergranular  cohe¬ 
sion:  the  strongest  interaction  is  exerted  between 
two  W3  atoms  across  the  GB  via  the  impurity 
atom  (see  Fig.  3).  We  shall  return  to  this  question 
in  the  following  section. 

The  GB-relaxation  calculations  also  show  that, 
as  in  the  case  of  the  Fe  GB,-^  the  interplanar 
separations  oscillate  as  a  function  of  distance  from 
the  GB,  the  deformation  waves  decaying  by  the 
10th- 12th  plane  away  from  the  GB.  An  interest¬ 
ing  feature  of  the  CL  GB  relaxation  is  that  the  dis¬ 
tance  between  the  2nd  and  3rd  planes  is  a  little 
over  a  half  of  the  (111)  interplanar  distance  in 
bulk  BCC  W  (0-550  A  versus  0-914  A).  Though 
W  does  not  Undergo  transformation  into  the  co¬ 
phase,  the  ‘misbalance’  in  interatomic  interaction 
arising  as  a  result  of  the  GB  in  an  attempt  to 
decrease  the  void  volume  between  two  W3  atoms 
across  the  GB  (see  Fig.  3)  results  in  the  tendency 
for  plane  3  nearly  to  collapse  onto  plane  2  (the 
&)-phase  configuration). 

The  impurity  atoms,  B,  C,  N,  and  O,  result  in 
some  ‘damping’  of  the  relaxation-deformation 
waves,  i.e.  decreasing  the  oscillation  amplitudes. 
This  damping  is  not  pronounced  for  B  and  C. 
Although  the  distance  between  planes  1  and  2 
(which  is  half  the  distance  between  two  atoms  in 
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planes  2  across  the  GB)  monotonically  increases, 
the  tendency  of  plane  3  to  coUapse  into  plane  2 
disappears:  in  the  progression  B  through  O,  the 
distance  between  planes  2  and  3  is  almost  equal  to 
that  in  the  bulk.  Except  for  this  distance,  the 
amplitude  of  the  deformation  waves  increases 
with  Al,  Si,  P,  and  S. 

The  analysis  of  GB  energetics  reveals  only  one 
aspect  of  an  impurity-embrittling  effect.  From  a 
thermodynamic  point  of  view,^‘^  the  impurity’s 
embrittling  potency  depends  on  the  difference 
between  the  free  energies  of  the  impurity’s  segre¬ 
gation  on  the  initial  GB  and  on  the  two  free 
surfaces  emerging  upon  fracture.  The  higher  the 
difference,  the  stronger  is  the  embrittling  potency 
of  the  impurity.  Thus,  in  order  to  have  a  complete 
thermodynamic  picture  of  the  embrittlement,  one 
should  be  able  to  analyze  the  energetics  of  the 
same  impurity  on  the  free  surface.  As  a  much  less 
rigorous  but  simpler  criterion,  the  sublimation- 
energy  differences  between  the  host  and  impurity 
atoms  were  calculated  by  Seah^  in  an  ideal-solu¬ 
tion  model  for  over  60  elements.  According  to 
Seah,"'  among  the  elements  analyzed,  only  B,  C, 
and  Os  may  be  cohesion  enhancers  in  W.  In  our 
more  rigorous  approach,  the  effects  of  impurities 
on  GB  stability  can  be  analyzed  by  simply  com¬ 
paring  the  GB-energy  differences,  AE,  between 
the  GB  with  impurities  and  the  CL  GB. 

The  corresponding  values  for  the  impurities 
discussed  are  plotted  in  Fig.  4.  One  can  see  that 
the  GB  stability  decreases  from  B  towards  O,  and 
the  energy  difference  becomes  positive  for  P  and 
S— the  strongest  embrittlers.  The  latter  means  that 
GBs  with  P  and  S  are  unstable  at  0  K. 


FIRST-PRINCIPLES  GB-SUPERCELL 
CALCULATIONS 

In  order  to  study  the  electronic  structure  of  the 
GB  (both  CL  and  that  with  an  impurity  atom  in  it) 
we  performed  a  series  of  LMTO  supercell  calcula¬ 
tions.  A  20-atom  supercell  was  used  as  a  model  of 
the  GB  (Fig.  5).  Again  as  in  our  semi-empirical 
calculations,  the  GB  is  modeled  by  the  succession 
of  ( 1 1 1 )  hexagonal  planes: 

CACBACBACBABCABCABCAC 

The  two  stacking  faults: . . .  C AC . . .  and  . . .  BAB . . . , 
imitate  two  (111)23  tilt  GBs.  The  filled  circles  in 
Fig.  5  show  impurity  atoms  (or  empty  spheres, 
if  the  GB  is  CL).  There  are  eight  planes  of  W 
atoms  between  the  two  A  GB  planes.  The 
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Fig.  5.  (a)  GB— schematic;  (b)  20-atom  hexagonal  supercell;  only  parts  adjacent  to  the  GBs  are  shown.  As  in  Fig.  1.  numbers 

designate  atoms  W1  through  W4;  (c)  view  along  [1 12]  direction  of  the  periodic  ( 1 1 1 )  plane  array;  small  and  large  circles  are  W 
atoms  in  alternating  (110)  planes.  Filled  circles— impurity  atoms  (or  a  vacancy  in  the  case  of  CL  GB). 


* 


immediate  impurity  environment  is  again  the  tri¬ 
gonal  prism  of  W  atoms  shown  in  Fig.  1. 

In  order  to  make  the  GB  model  more  realistic, 
the  interplanar  distances  were  taken  to  be  equal  to 
those  obtained  from  the  above  semi-empirical 
relaxation  calculations  for  planes  2  through  5  (on 
both  sides  of  the  GB  planes);  the  distance  between 
the  two  equivalent  planes  in  the  middle  of  the 
supercell  was  set  equal  to  the  interplanar  distance 
in  the  bulk  (0-914  A). 


When  a  metalloid  atom  is  added  into  a  transi¬ 
tion-metal  crystal  lattice,  two  effects  are  pro¬ 
duced.-^  First,  a  covalent  bond  between  the 
impurity  atom  and  the  host  transition  atom  is 
formed,  and  then  the  crystal  lattice  is  expanded 
and  deformed  (in  fact,  both  processes  occur 
simultaneously).  A  similar  situation  takes  place  in 
a  W  GB. 

Figure  6  shows  the  site-projected  electronic 
densities  of  states  (DOSs)  for  atoms  W1  and  W3 


Fig.  6.  The  site-projected  DOSs  (states/Ry)  for  atoms  W1  and  W3.  Zero  energy  here  and  in  Fig.  8  corresponds  to  the  £p  of 

bulk  BCC  W. 
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for  the  CL  GB  and  GBs  with  different  impurities. 
The  lower  (‘negative')  parts  of  the  plots  show  the 
site-projected  DOSs  within  the  atomic  spheres  of 
the  corresponding  impurities;  they  are  identical  in 
plots  for  W1  and  W3  atoms.  The  plots  in  Fig.  6 
actually  represent  the  GB  electronic  band  struc¬ 
ture,  and  allow  an  analysis  of  interatomic  bonding. 

When  an  impurity  atom  is  immersed  into  the 
electron-atom  system  of  the  host,  its  electrons 
become  part  of  the  whole  system.  The  impurity’s 
electronic  states  hybridize  with  the  electronic 
bands  of  the  host  metal,  which  results  in  the  for¬ 
mation  of  covalent  bonds  between  the  impurity 
and  the  host  atoms.  In  fact,  the  whole  electron- 
charge-density  distribution  becomes  affected,  and 
interactions  among  all  the  atoms  are  disturbed. 

As  is  typically  the  case  for  metalloid  im¬ 
purities,'"’  the  impurity  A-electrons  (25  for  B 
through  O,  and  35  for  P  and  S)  form  in  W  the  nar¬ 
row  impurity  bands  below  the  bottom  of  the 
metal-valence  band.  Even  though  the  correspond¬ 
ing  5-levels  in  the  free  atoms  lie  well  below  the  W 
p-  and  <3f-bands,  upon  alloying,  the  5  electrons  do 
hybridize  with  those  bands,  resulting  in  strongly 
localized  but  rather  weak  ionic-type  bonds. 

The  impurity’s  p-electrons  with  the  energy 
right  inside  the  W  (/-bands,  create  pronounced 
covalent  bonds.  The  impurity-bonding  states  of 
predominantly  p-type  hybridize  with  W  mostly  d- 
and  also  p-states.  Comparison  of  the  DOS  plots 
for  W 1  and  W3  atoms  show  that  the  •-W 1  hybri¬ 
dization  in  the  GB  plane  is  stronger  than  #-W3 
across  the  GB,  in  spite  of  the  fact  that  the  former 
distances  are  larger  than  the  latter  (Table  2).  This 
means  that  the  W  atoms  in  the  GB  plane  are  more 
strongly  bound  than  those  across  the  GB.  As  one 
can  see,  the  impurity  s-p  bands  are  shifting  to¬ 
wards  the  W  (/-band  bottom,  thus  ‘switching  off’ 
the  W  bonding  (/-electrons  in  the  upper  part  of 
the  (/-bands.  As  a  result,  the  hybridization 
becomes  weaker  in  the  row  B  through  O.  In  the 
case  of  B,  a  significant  part  of  the  W  (/-states  is 
involved.  In  fact,  one  can  even  speculate  that  the 
•-W3  hybridization  is  the  strongest,  since  a  pro¬ 
nounced  peak  of  W  (/-states  is  involved.  For  O,  P, 
and  S,  the  hybridization  of  electronic  (/-states  at 
the  bottom  of  the  (/-band  (the  peaks  around 
—  600  mRy)  within  the  atom  W3  sphere  almost 
disappears,  while  it  is  still  moderate  within  the 
atom  W 1  spheres. 

Another  important  trend  can  also  be  seen. 
Beginning  with  N  and  beyond,  some  anti-bonding 
states  (in  both  W  and  •)— to  the  right  of  the 
‘troughs’  around  the  Fermi  energy,  £"p— are  pro¬ 


gressively  filled,  which  results  in  weakening  of  the 
interatomic  bonds.  In  terms  of  the  charge  density, 
‘bonding’  means  a  pile-up  of  electronic  charge  in 
the  space  between  the  atoms,  whereas  the  ‘anti- 
bonding’  states  result  in  decreasing  the  charge 
density.  Figure  7  compares  the  electron-charge 
densities  in  atomic  spheres  of  B  and  O.  The 
arrows  show  the  density  values  corresponding  to 
the  so-called  ‘muffin-tin’  radii,  the  radii  Of  touch¬ 
ing  hard  spheres  representing  a  B  or  O  and  W3 
atom  (the  nearest  neighbor  to  the  impurity  atoms). 
One  can  see  that  the  electron  density  of  B  at  the 
‘point  of  contact’  with  atom  W3  is  1-5  times  as 
great  as  that  of  O,  which  thus  suggests  a  higher 
‘bonding  capacity’  of  B,  as  compared  with  that  of 
O. 

Finally,  it  is  interesting  to  compare  the  total 
DOSs  of  a  CL  GB,  GBs  with  impurities  and  bulk 
BCC  W.  Figure  8  shows  the  plots.  As  in  the  site- 
projected  DOSs,  the  GBs  with  impurities  demon¬ 
strate  strongly  localized  impurity  bands.  Each  of 
them  contain  exactly  two  electrons,  as  should  be 
the  case,  since  the  corresponding  Brillouin  zones 
are  completely  filled.  As  was  mentioned  earlier, 
the  anti-bonding  states  (to  the  right  of  the  troughs 
around  100  mRy),  being  almost  completely 
unfilled  in  bulk  W,  begin  progressively  to  fill  in  the 
GBs.  The  Fermi  energy,  Ep,  sits  on  peaks,  their 
heights  increasing  from  B  through  S.  In  fact,  the 
Ep  peaks  in  the  CL  GB  is  higher  than  that  for  B.  A 
relative  DOS  value  at  Ep  is  known  to  be  an  indica¬ 
tion  of  the  system  stability.  The  higher  the  value, 
the  lower  is  the  stability.  From  this  point  of  view, 
the  stability  of  the  GB  is  lower  than  that  of  the 
bulk  W,  whereas  B  improves  the  stability  as  com¬ 
pared  with  that  of  the  CL  GB. 
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Fig.  7.  Electron-charge  densities  in  atomic  spheres  of  B 
and  O. 
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Fig.  8.  Total  DOSs  for  BCC  W,  CL  GB,  and  GB  with 
impurities.  Arrows  show  the  Fermi  energies. 


CONCLUSIONS 

Both  the  semi-empirical  and  first-principles 
analyses  have  shown  that  B  in  the  W  GB  plays  a 
dual  role.  First  of  all,  owing  to  the  site-competi¬ 
tion  effect,  it  tends  to  displace  the  other  impurity 
atoms  off  the  GB,  which  thus  ‘cleanses’  it.  At  the 
same  time,  B  enhances  the  intergranular  cohesion, 
which  thus  improves  resistance  to  brittle  fracture. 
A  simple  estimate  shows  that  10-50  ppm  of 
impurity  atoms  will  saturate  the  GBs  in  W.  Ideally, 
the  same  amount  of  B  would  be  sufficient  to 
improve  the  ductility  significantly.  However,  the 
above  analysis  disregards  a  possible  chemical 
activity  of  B,  e.g.  forming  boron  oxides  or  tung¬ 
sten  borides.  The  latter  would  require  introducing 
a  multiplying  factor  to  correct  for  B  ‘consumed’  as 
boron  compounds.  Alloying  W  with  B  in  quan¬ 
tities  10“  15  times  as  great  did  result  in  a  sig¬ 
nificant  ( 150®C)  drop  in  DBTT.^’  This  effect  was 
attributed  to  gettering  O  by  forming  boron  oxides. 
In  any  case,  microalloying  W  with  B  is  extremely 
promising.  The  experimental  work  directed  at 


elucidating  the  various  aspects  of  microalloying  W 
is  currently  in  progress  at  the  US  Army  Research 
Laboratory. 
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